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ABSTRACT 

Coronal elemental abundances are known to deviate from the photospheric values of their parent star, with the degree of deviation 
depending on the first ionization potential (FIP). This study focuses on the coronal composition of stars with supersolar photospheric 
abundances. We present the coronal abundances of six such stars: 11 LMi, l Hor, HR 7291, r Boo, and a Cen A and B. These stars all 
have high-statistics X-ray spectra, three of which are presented for the first time. The abundances we measured were obtained using 
the line-resolved spectra of the Reflection Grating Spectrometer (RGS) in conjunction with the higher throughput EPIC-pn camera 
spectra onboard the XMM-Newton observatory. A collisionally ionized plasma model with two or three temperature components is 
found to represent the spectra well. All elements are found to be consistently depleted in the coronae compared to their respective 
photospheres. For 11 LMi and r Boo no FIP effect is present, while l Hor, HR 7291, and a Cen A and B show a clear FIP trend. These 
conclusions hold whether the comparison is made with solar abundances or the individual stellar abundances. Unlike the solar corona, 
where low-FIP elements are enriched, in these stars the FIP effect is consistently due to a depletion of high-FIP elements with respect 
to actual photospheric abundances. A comparison with solar (instead of stellar) abundances yields the same fractionation trend as on 
the Sun. In both cases, a similar FIP bias is inferred, but different fractionation mechanisms need to be invoked. 
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. Introduction 


The active coronae of cool stars can be up to five orders of mag¬ 
nitude more luminous and an order of magnitude hotter than our 
own solar corona. Like the Sun, they emit both extreme UV 
(EUV) and soft X-rays. X-ray spectroscopy allows for partic¬ 
ularly deep insights into the hot thermal content and elemen¬ 
tal abundances of the corona. For comprehensive reviews see 
Giidel & Naze| ( |2009| l and |Testa| ( |2010| >. Solar coronal abun¬ 
dances are different from the solar photospheric composition. 
The abundance patterns depend on the first ionization potential 
(FIP) ( Meyer| 1 985[ Feldman|1992| ). Elements with FIP<10 eV 
(low FIP) are enriched in the corona relative to elements with 
FIP>10 eV (high FIP). Averaging over the entire solar disk, the 
enrichment factor is about 4 ( Feldman & Widing|2002 and ref¬ 
erences therein). Currently, there is no generally accepted model 
that explains the FIP bias (but see Laming 2004| |Telleschi et al. 


2005 1 Laming|2012| l. Nevertheless, the observed reality of abun¬ 
dance fractionation is undisputed, although the absolute normal¬ 
ization of the solar coronal abundances is still being debated (cf., 
Schmelz et al.|2012] ). 


The FIP effect became even more puzzling when high- 
resolution X-ray spectra from XMM-Newton and Chandra revealed 
that active stellar coronae do not follow the solar FIP pattern. 
In some cases, the high-FIP elements are even enriched com¬ 
pared to the low-FIP ones, an effect labeled inverse-FIP (IFIP) 
effect ( |Brinkman et al.|2001| >. Later studies revealed that FIP and 
IFIP biases are correlated with coronal activity (e.g.. X-ray lu¬ 
minosity) and age: Highly active stars show an IFIP effect, while 
less active coronae have a solar FIP bias (Aud ard et al.| [2003; 
Telleschi et al. 2005), implying a transition on stellar evolution 
time-scales of Gyr. Intervening intermediate-activity stars seem 
to exhibit a relatively flat dependence on FIP. [Ball et al.|(|2005) 


even suggested a U-shaped dependence, although such subtleties 
depend on one or two elements at most. Some more recent ob¬ 
servations of stellar coronae have found that there are deviations 
from this behavior for some stars, however, and that there could 
be additional factors that determine coronal abundances, for ex¬ 
ample, spectral type ( Giidel et al.||2007| Wood & Linsky|2010[ l. 
In particular, the influence of stellar metallicity on stellar coro¬ 
nal properties is still unclear, since most of the stars with well- 
studied X-ray spectra have photospheric abundances similar (or 
assumed to be) to the solar photospheric values. 

Moreover, as is the case on the Sun, different parts of a sin¬ 
gle corona may exhibit different abundances (Sanz-Forcada et al. 
2004| [Nordon et aL 2013). Wood & Linsky (2006) found that 
even though two stars may have very similar properties and ac¬ 
tivity levels (and being visual binary companions), 70 Oph A and 
70 Oph B did not exhibit the same coronal abundance patterns. 


Given the importance of the heavy elements in the cooling 
of plasmas with temperatures in the range from 10 6 to 10 75 K, 
it is expected that metallicity could play a significant role in de¬ 
termining the thermal structure (the emission measure distribu¬ 
tion, or EMD) of the coronae of such stars. For example, one 
might expect that in the absence of fractionization mechanisms, 
coronae with enhanced (depleted) heavy elements would be sig¬ 
nificantly cooler (hotter) than coronae with solar abundances. It 
is thus of considerable interest to study the coronal spectra of 
stars whose photospheric abundances differ the most from so¬ 
lar, so as to study how varying the source plasma abundances 
affects the properties of the coronal plasma that ultimately must 
originate from the underlying photosphere. In the present anal¬ 
ysis we select several stars that are metal enriched with respect 
to solar from a mean factor of 1.5 to 2.7. We seek to determine 
whether coronal abundances behave in a FIP or IFIP behavior. 
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and whether this behavior is a result of depletion or enrichment 
in either low-FIP or high-FIP elements. 

The present study aims to test whether the FIP (or IFIP) ef¬ 
fect can be accurately determined relative to stellar abundances. 
While the standard assumption of solar photospheric abundances 
may be valid in many cases, a true FIP trend requires the knowl¬ 
edge of true photospheric abundances, which was stressed by 
Sanz-Forcada et al. (2004). Photospheric abundances are diffi¬ 
cult to measure and to some extent are model dependent, while 
coronal abundances are comparatively easily measured by count¬ 
ing X-ray photons in emission lines. The result is that a surpris¬ 
ingly small number of bright coronal sources also have reliable 
photospheric abundances. We proposed to observe 11 LMi after 
searching the literature for bright X-ray sources that have reli¬ 
able photospheric abundances that are markedly different than 
solar. High-precision photospheric abundances of stars bright 
enough in X-ray such that well-exposed high-dispersion X-ray 
spectra can be obtained are few and far between in the literature. 

Since different references for solar abundances are used by 
different authors, it is important to normalize all measurements 
to a standard solar reference. Abundances presented in this paper 
are normalized to the solar reference of Asplund et al. (2009 1 . 


minosities span two orders of magnitude, from 1 ()-' to 10 29 erg 
s' 1 . Although different photospheric abundance measurements 
for these object exist in the literature, they are predominantly 
supersolar. 

The coronal abundances of 11 Leonis Minoris (11 LMi) , 
z Horologii (z Hor), and HR 7291 are presented here for the first 
time, while measurements for the other three, t Bootis (r Boo 
) as well as a Centauri (a Cen) A and B, are already available 
in the literature. The coronal abundances of t Boo and a Cen 
A, B were extensively analyzed by Maggio et al. ( 201 1|) and by 
Raassen et al.| ( 2003]> respectively. 


Maggio et al. (|2011 1 found 


( 119911 , with an X-ray luminosity much lower than that of the RS 
CVn class of close binaries that usually exhibit Lx = 10’° 52 erg 
s -1 . It is classified in the spectral class of G8+V in line with an 
G8 IV-V object. There is a much fainter (AV = 7.7 mag) visual 
binary companion 5" away from the G star, probably a mid-M 
dwarf with a binary orbit of 201 years (see jMalkov et al|2012 1 . 
This companion is not the X-ray source since the ROSAT HRl 
source lies within l"of the position of the G star. Soubiran & Gi¬ 
rard (|2005j> considered 11 LMi to probably be a (95%) thin-disk 


star, consistent with its enhanced rotation rate and activity level 
compared to the Sun. 

2.2.2. z Hor 

z Hor (HD 17051) is a planet-hosting star that is dim in X-rays. 


The orbit of the exoplanet is approximately 311 days (Sanz- 
|Forcada et al. 2013|>. In addition, it has a 1.6 year magnetic activ¬ 
ity cycle ([Metcalfe et al.|2010|l. X-ray observations are detailed 


Sanz-Forcada et al. (2013) where a XMM-Newton campaign 


was used to observe the activity cycle. 


2.2.3. HR 7291 


2. Sample 

2.1. Object selection 

Objects for this study were chosen according to two criteria - 
high photospheric metallicity, and high-statistics X-ray spectra. 
Sources were filtered from Gonzalez & Laws (2007) and A llende| 
Prieto et al. (2004). First, objects with supersolar abundances of 
C, O, Mg, Si, Fe, and Ni were selected if at least four of these 
sustained Az/H > 1.1. From these, all objects with X-ray grating 
observations of at least 1000 photon counts were selected, leav¬ 
ing five stars, two of which are in a binary system. In addition, 
11 LMi was observed specifically in preparation for this analysis 
(see jDrake et al.|2010[ ). A summary of objects, their references, 
and some of their properties is given in Table [I] 

The sample includes six nearby ( <27 pc ) stars with super¬ 
solar metallicities ranging from approximately 1.5 to 2.5 solar. 
Spectral types of the sample are between K1 and G8. X-ray lu- 


HR 7291 is an F8V star with a planet with an orbital period of 3 
days. X-ray observations in this case are detailed in Scandariato 
et al. (|2013|, as part of a campaign to search for a planet-induced 


activity cycle. 


2.2.4. r Boo 


a flat FIP dependence, while Raassen et al. (2003) measured a 
solar-like FIP behavior. 


2.2. Additional notes on objects 
2.2.1. 11 LMi 

A nearby moderately X-ray active star, 11 LMi has been erro¬ 
neously classified as an RS CVn binary by Simbad - probably 
on the basis of its detection as a rotationally modulated (18.0 
days) variable star ( Skiff & Lockwood||1986|). In fact, 11 LMi 
is not a spectroscopic binary according to Duquennoy & Mayor 


This X-ray bright star hosts a planet with a 3.3-day period. [Mag-| 
|gio et al!] ( |201 1 > analyzed the FIP behavior of this object. No 
clear FIP (or IFIP) trend was observed. Instead, a flat depen¬ 
dence relative to both solar and actual photospheric abundances 
was found. In addition, there is a factor 4 of depletion across all 
elements observed in the corona relative to the photosphere. 

2.2.5. cr Cen binary 

This spectroscopic binary is part of the nearest stellar system to 
the Sun (1.34 pc). Both stars are bright X-ray sources. Raassen 
|et al.| ( [2003l > observed a FIP effect for both stars, relative to solar 
abundances. 


2.3. Observations 

The log of observations for the present sample is listed in Ta¬ 
ble [2] Count rates were averaged and exposure durations were 
summed. Four of these stars were observed in X-rays by XMM- 
Newton , while a Cen A and B were observed using Chandra. 

The observations of 11 LMi , r Boo , a Cen A, and B have 
sufficiently long exposures that result in high-quality grating 
spectra (see Maggio et al.|201 1 |Drake et al.|2010; Raassen et al. 


|2003| for details), making for reliable abundance measurements. 
On the other hand, z Hor and HR7291 were both observed inter¬ 
mittently during long observation sequences. As such, the pho¬ 
ton count in each individual observation is low (around 200), and 
many observations exist (Sanz-Forcad a et al.|2013| [Sc andariato 
et al.||2013). This makes the high-resolution spectra of the Re¬ 


flection Grating Spectrometer (RGS) noisy, and it is difficult to 
determine lines. Thus for these two objects the CCD-based Eu¬ 
ropean Photon Imaging Camera (EPIC) is more reliable for an 
abundance analysis (see Sect.|3.2.2|i. 
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Table 1: Object details 


Object 

Spectral Type 

Distance 

pc 

Metallicipri 

AJA q 

Reference^] 

Total Grating Count^] 

L x 

10 28 ergs~' 

11 FMi 

G8IV-V 

11.4 

2.06 

1,2,3,4,5 

6588 

3.2 

l Hor 

GOV 

17 

1.47 

6,7,8 

7785 

1.6 

HR 7291 

F8V 

27 

1.49 

6,7 

1345 

3.7 

r Boo 

F7V 

15.6 

2.37 

6,9 

9741 

10 

a Cen A 

K1V 

1.34 

2.67 

10,11 

11840 

0.1 

a Cen B 

GOV 

1.34 

2.15 

10,11 

11360 

0.1 


1 Average abundance in photosphere (averaged over all references), with respect to solar 

Asplund et al. 

2009 





2 1 

Soubiran & Girard 

(20051, 2 

Pmgniel 

(20111,3 

Ramirez 

2007 

i,4 

^uck & Heiter 

20 

06 

, 5 Milone 

2011 

1,6 

Gonzalez & Laws 

2007 


7 

Biazzo et al. i 

2012 

8 

Bond et al. i 20( 

)6 ', 9 Takeda et al. 

2001 

ic 

Allende Prieto et af 

(2004), 11 

Porto 

de 

Vlello et al. 

(2008) 



3 Zeroth order for Chandra 


3. Coronal abundance measurements 

We measured the coronal abundances for three of the targets: 
11 LMi , l Hor, and HR 7291. These targets do not have pub¬ 
lished detailed X-ray abundances. 


3.1. Methodology 


We used the collisional ionization equilibrium (CIE) code VApec 
( Smith et al. 2001 1 within the standard XSPEC1 10 software 
package. Models were selected according to the best-fitting re¬ 
duced x 2 goodness test, and when several models were consid¬ 
ered, the F-Test was used to determine the statistical robustness 
of the model selection. All fitting results, reduced^ 2 , F-test, and 
uncertainty values were computed using XSPEC. All spectral 
plots in the paper show data rebinned to a minimal significance 
of 3cr for each data point, but not more than three bins, so as to 
retain the high resolution in the figures. However, the spectral 
fits use the original non-binned higher resolution. 


In general, the RGS spectra would be preferred because of 
the high spectral resolution and the ability to visually inspect 
emission lines. However, meaningful results can only be ob¬ 
tained for 11 FMi (see Sect. |3.2.2[ i. When no constraint could 
be obtained using the RGS, the abundances were measured us¬ 
ing the EPIC-pn spectra. EPIC-pn is also advantageous for ob¬ 
serving Mg and Si, which are visibly discernible and are gener¬ 
ally lacking in the RGS spectra. We concentrated on the well- 
exposed EPIC-pn spectrum since it has more than enough counts, 
and hence we did not use data from the other EPIC-MOS detec¬ 
tors on XMM-Newton. For i Hor and HR 7291 none of the RGS 
data were found useful (again, see Sect. |3.2.2| ). 

Joint fits between RGS and EPIC-pn were considered but dis¬ 
carded. Since the EPIC-pn provides such a high photon count, 
the fit is dominated by the EPIC-pn alone and the high spectral 
resolution of the RGS is lost. Since both i Hor and HR 7291 
have multiple low-count observations, RGS data are highly noisy 
and only Fe/O ratios can be constrained. Therefore only EPIC-pn 
was used for these targets. A joint fit of all observations was 
attempted, but significant differences in flux levels between ob¬ 
servations rules out a single model that fits all data well (reduced 
X 2 > 2). Instead, each observation was fit separately and results 
were averaged (see Sect. |3. 2.2}. 


4 http://heasarc.nasa.gov/xanadu/xspec/ 


3.1.1. RGS analysis 


When modeling the RGS spectrum (11 FMi only), we used 
wavelengths between 7 - 37 A where the S/N is high enough to 
visibly detect lines. The free parameters of each thermal compo¬ 
nent are the temperature T and the elemental abundances A z of 
the plasma, as well as the emission measure (EM = f n e n H dV). 
Testing for multiple T components, for IT,2T,3T, and 47’ the 
reduced x 2 is minimal for the 37’ model with a value of 1.08. 
In addition, we conducted an F-test to check for statistical ad¬ 
vantage of adding components. As expected, the model with 3 T 
model is statistically advantageous (null hypothesis probability 
of < 10 ' when compared to any other model). We thus deem 
the 3 T model to be the most appropriate approximation of what 
is probably a continuous EMD. All models with two or more 
temperatures yield consistent relative abundances to within the 
errors, which increases our confidence in the derived values. 


Abundances were held constant between temperature com¬ 
ponents assuming abundance dispersion does not change along 
the EMD. Since no hydrogen lines are available in X-ray spectra, 
it is necessary to freeze the abundance of one of the observed ele¬ 
ments to avoid the inherent degeneracy between the overall metal 
abundances and the EM. This limits the measurement to relative 
abundances of the observed elements. We chose to fix the Fe 
abundance to its solar value (again, Asplund et al .|2009 1, as Fe- 
F lines originate from a wide range of T. Hence, we were able 
to fit for the abundances of C, N, O, Ne, and Ni while all other 
elements were set to their solar values. Indeed, the attempt to fit 
for the Mg and Si abundances failed, which forced us to com¬ 
pletely freeze their values to solar as well. Nevertheless, the 3 T 
collisional ionization model reproduces all featured RGS spec¬ 
tral lines well, providing robust results for our free abundance 
parameters and their uncertainties. Moreover, setting the Mg and 
Si values to arbitrarily high or low values has little to no effect on 
the fit, and no meaningful constraint on these abundances could 
be derived from the data. 


3.1.2. EPIC-pn 

Since the absolute coronal abundances (relative to H) are key to 
this study because of the supersolar abundances in the photo¬ 
sphere (Sect. [TJ and since the RGS falls short in providing ab¬ 
solute values, we must constrain absolute abundances with the 
EPIC-pn camera. While not as accurate as the RGS in determin¬ 
ing line positions, the high sensitivity of the EPIC-pn is advan¬ 
tageous for measuring the bremsstrahlung continuum, which is 
mostly due to H. We thus fit the EPIC-pn spectrum without fixing 
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any of the observable elements to check whether we could obtain 
absolute abundances that are not degenerate with the normaliza¬ 
tion (EM). Since the effective area curve of EPIC-pn covers the 
band of ~ 1.5 - 30 A (0.4 - 10 keV), which is somewhat differ¬ 
ent from that of the RGS, the model for fitting its spectrum may 
have a less detailed thermal structure. We tested both 2 T and 3 T 
models, as suggested by the RGS fits, (for i Hor and HR 7291 as 
well). The third T component does not improve the fit (reduced 
X 2 ), nor does it reduce the uncertainties on the model param¬ 
eter values. Consequently, we preferred the simpler 2 T model, 
and all EPIC-pn results are given for this model. The Mg and Si 
emission lines, which are usually lacking in most RGS spectra, 
are clearly seen around 9 A and 6.6 A in all EPIC-pn spectra we 
analyzed. 


3.2. Results 


3.2.1. 11 LMi 

The full RGS spectrum of 11 LMi is presented in Fig. [T] along 
with the residual plot beneath. An inspection of the spectrum 
clearly indicates the presence of the usual main coronal emission 
lines (e.g., Behar et al. 2001!, with the Fe L-shell dominating 
between 12 - 18 A, O K-shell at 19 - 23 A, and C K-shell at 
~ 34 A being strong as well. The Ne-K (12 - 14 A) and N-K 
(~ 25 A) lines are weaker. Still, they are easily seen. At lower 
wavelengths, Mg-K is just barely detected around 8 - 9 A, while 
Si-K (6 - 7 A) is not detected at all and is not included in Fig. Ill 
The best-fit parameters of the model are presented in Table B] 
along with their lcr uncertainties. The abundance ratios of Fe are 
fixed to solar in the RGS analysis, but in the EPIC-pn analysis we 
find they are consistent with solar values at the lcr level. 

The EPIC-pn spectrum and best-fit model are presented in 
Fig. 2 and are included in Table [3] For the most part, all values 
are consistent between RGS and EPIC-pn to within the lcr level 
(except Ni/H at 2cr), but clearly, some abundance measurements 
are not very constraining. Certainly, the low spectral resolution 
available with EPIC-pn and the low S/N available with the RGS 
could be blamed for this. We used the RGS absolute abundances 
of C, N, O, Ne, and Ni, and the EPIC-pn for Mg, Si, and Fe. 

We next wish to ensure that the best-fit solution is not de¬ 
generate between its overall metal abundances (metallicity) and 
the total EM. To obtain a confidence contour for the metallicity 
in the corona of 11 LMi, we tied all best-fit abundance ratios to 
Fe and also tied the EM ratio of the two components (temper¬ 
atures were fixed), effectively leaving only two free parameters 
- namely the metallicity and EM. The resulting confidence con¬ 
tours are presented in Fig. [3] The Fe abundance (with all metals 
tied to their EPIC-pn fitted ratio to Fe) is tightly constrained be¬ 
tween 0.8 - 1.1 to 3cr. This result gives a tighter constraint, yet 
is fully consistent with the 0.8 - 1.5 (lcr) Fe abundance of the 
full model, where all abundances were allowed to vary. 


3.2.2. i Hor and HR 7291 

For l Hor and HR 7291 only the EPIC-pn spectra were used. Tem¬ 
perature variability between observations was the goal and is in¬ 
deed present during these campaigns. Consequently, a simulta¬ 
neous fit to all observations with a single 2T CIE model yields an 
unacceptable fit (reduced^ 2 > 3). Therefore, each of the 14 (15) 
observations of i Hor (HR 7291) were fitted separately. Since we 
do not find significant variations of the abundances to within 2cr 
and the vast majority are consistent to lcr, we averaged the best- 
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Fig. 1: RGS spectrum of 11 LMi between 7 - 37 A (data points) 
and best-fit model (red solid line). Residuals are plotted in the 
bottom panel. The data in the figure are rebinned to 3cr, but not 
more than three bins. 


fit abundance values across observations. Admittedly, a few ob¬ 
servations are clearly of poor quality, as evident by their reduced 
X 1 > 2 fits or their very large parameter uncertainties. Hence, we 
cannot rule out abundance variability. Nevertheless, to avoid out¬ 
liers due to low statistics, the few observations with abundances 
that are more than lcr different from the average were discarded 
(two for c Hor, four for HR 7291). If the abundances are truly 
constant, we significantly improve the constraints on the fitted 
abundances by averaging the results. The standard error was cal¬ 
culated for each elemental abundance and is consistently lower 
than the statistical errors of each fit (by a factor of 2 at most). Of 
course, this result is only valid if CIE is prevalent in the plasma 
for each individual observation. 

The results are presented in Table[4]for the best-fitting mod¬ 
els with average reduced^ 2 = 1.0 for both l Hor and HR 7291, 
including observations with both high and low reduced^ 2 (>1.2 
and <1). Clearly, abundances averaged from all fits provide well- 
constrained values for all elements, although some could not be 
resolved by EPIC-pn or are missing completely, for example, N. 
The same method as described in Sect. l3.2.H was used to obtain 
contour plots between metallicity and EM. We present here a 
single contour plot for example, since most enclose similar val¬ 
ues in parameter space and our final constraints are much bet¬ 
ter than that of a single observation. Spectra and contour plots 
are presented in Figs. [4] and [5] The fact that the Fe abundance 
(representing metallicity) is tightly constrained to 3<x in all three 
analyzed objects to a range of 0.3 - 0.8 in solar units boosts our 
confidence in the absolute abundances with respect to H of all 
elements. 

4. Coronal abundances from the literature 

Coronal abundances of r Boo and a Cen A and B are only quoted 
here and not measured again, since high-quality measurements 
are already available. In Table[5]we list these and the above mea¬ 
sured abundances and compare them with their respective photo- 
spheric abundances taken from the literature (see Table[Ii. Since 
N and Ne are only measured in the coronae and not in the pho¬ 
tospheres, an effective comparison can only be made for five el¬ 
ements - C, O, Mg, Si, and Fe. All abundances are given relative 
to Asplund et al. (|2009[>. Note that photospheric ranges repre- 
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Table 2: Observation log 


Object 

Dattj^j 

Instruments 

Telescope 

Total Duration 
ks 

11 LMi 

09/5/6-7 

RGS, EPIC-pn 

XMM-Newton 

85 

i Hor 

11/5/16-13/2/3 

RGS, EPIC-pn 

XMM-Newton 

77 

HR 7291 

09/9/21-09/10/27 

RGS, EPIC-pn 

XMM-Newton 

24 

r Boo 

03/6/24 

RGS, EPIC-pn 

XMM-Newton 

65 

a Cen A 

99/12/25 

LETG-HRC-S 

Chandra 

82 

a Cen B 

99/12/25 

LETG-HRC-S 

Chandra 

82 


1 Date range represents general times of observations, not one continuous observation. 


Table 3:11 LMi summary of fitted parameters 


Parameter 

EPIC-pn Value 

Uncertainty 

lcr 

RGS Value 

Uncertainty 

lcr 

x 2 

273 


5375 


DOF 

265 


5002 


kTi [keV] 

0.23 

0.22-0.24 

0.14 

0.11-0.19 

kT 2 [keV] 

0.57 

0.55-0.60 

0.38 

0.34-0.41 

kT 3 [keV] 



0.69 

0.57-0.83 

C [solar] 

0.1 

0.0-4.2 

1.1 

0.6-1.9 

N [solar] 

1.5 

1.0-3.8 

0.6 

0.2-1.2 

O [solar] 

0.4 

0.3-0.8 

0.8 

0.6-1.0 

Ne [solar] 

0.9 

0.5-2.0 

0.8 

0.6-1.0 

Mg [solar] 

0.9 

0.7-1.5 



Si [solar] 

0.6 

0.4-0.9 



Fe [solar] 

1.0 

0.8-1.5 

1.0 

fixed 

Ni [solar] 

4.4 

2.0-7.1 

0.8 

0.0-1.9 


Table 4: t Hor and HR 7291 best-fit parameters 


Parameter 

Value 

l Hor 

Uncertainty 

lcr 

HR7291 

Value Uncertainty 
lcr 

Avg.* 2 

122 


117 


Avg. DOF 

116 


118 


kT x [keV] 

0.26 

0.23-0.30 



kT 2 [keV] 

0.75 

0.68-0.81 



O [solar] 

0.4 

0.3-0.5 

0.25 

0.2-0.3 

Ne [solar] 

1.1 

0.9-1.3 

0.4 

0.3-0.5 

Mg [solarl 

1.9 

1.1-22 

0.9 

0.8-1.0 

Si [solar] 

1.5 

1.2-1.7 

0.3 

0.2-0.4 

Fe [solar] 

1.3 

1.2-1.4 

0.8 

0.7-0.9 


sent those values found in the literature, whereas coronal ranges 
are the measured lcr uncertainties. Overall, coronal abundances 
are significantly lower than their supersolar photospheric coun¬ 
terparts. In the next section the FIP trends of all these stars are 
discussed. 


5. Discussion 

5.1. FIP effect 

The dependence of abundances on FIP is plotted in Fig. [7] Abun¬ 
dances are normalized to solar (Asplund et al. |2009j ), and photo- 
spheric values are shown for comparison. The FIP behavior may 
thus be compared to both solar and actual abundances. While the 
uncertainties of the coronal abundances in the figure are those 


measured in the present analysis or referenced in Sect. 2.1 the 


uncertainties of the photospheric data only reflect the range of 
values in the literature. The lack of real errors on the photo¬ 
spheric values somewhat hinders a meaningful comparison, al¬ 
though by combining several literary sources we obtain a sense 
of the constraint on these abundances. 

The first obvious result observed is that depletion of elements 
differs significantly when comparing to solar abundances instead 
of using actual photospheric abundances. Furthermore, high-FIP 
elements are strongly depleted in all stars in the sample. Thus, 
objects displaying a FIP effect show little to no depletion of 
low-FIP elements (unlike the solar corona, which displays an 
increase of low-FIP elements), and objects displaying no FIP ef¬ 
fect show a more even depletion (11 LMi, r Boo ). For a Cen, a 
comparison with solar and actual photospheric abundances leads 
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Wavelength (A) 


Fig. 2: EPIC-pn spectrum of 11 LMi with best-fit model (in red). 
Mg and Si emission lines are seen around 9 A and 6.6 A. The 
high statistics allows us to obtain absolute abundances by con¬ 
straining the line emission (metals) with respect to the contin¬ 
uum (mostly H). Residuals are plotted in the bottom panel. The 
data in the figures are rebinned to 3cr, but not more than three 
bins. 


O 



Fig. 3: Contours for 11 LMi of overall normalization represent¬ 
ing EM (ratio of two thermal components tied) and Fe abun¬ 
dance, while ratios of all other metal abundances to Fe (Az/Ap e ) 
are set to their EPIC-pn best-fit values. The ability to obtain ab¬ 
solute metal abundances with confidence to the level of 3cr is 
demonstrated. 


to the same FIP behavior, but to differing conclusions regarding 
the origin of the effect. 

11 LMi shows no FIP dependence with respect to solar abun¬ 
dances. When comparing with the supersolar photospheric abun¬ 
dances, most significantly Si, but also Ni and possibly Mg are de¬ 
pleted, which may give the impression of an inverse FIP (IFIP) 
effect. However, Fe, which is a low-FIP element, is not signifi¬ 
cantly depleted, while O, a typical high-FIP element, is depleted, 
both of which argue against an IFIP effect. 

In i Hor the low-FIP elements do not significantly deplete 
with respect to the photosphere, although they appear to be su- 
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Fig. 4: l Hor (top) and HR 7921 (bottom) EPIC-pn spectra, model, 
and residuals. 


persolar. O in the corona is significantly depleted with respect 
to the supersolar photospheric value, producing a solar-like FIP 
effect. 

HR 7291 shows a rather flat depletion across the FIP axis. 
There seems to be an enhanced depletion with increasing FIP 
that creates a solar-like FIP trend as in i Hor . In this case, the 
depletion of Ne reinforces this effect, r Boo is similar to 11 LMi 
as all elements are depleted relative to the photosphere with no 
I/FIP effect evident. Depletion in this case is much stronger than 
in the actual photospheric abundances. 

The two a Cen stars exhibit a stronger FIP effect in their 
coronae when compared relative to their actual photospheric 
abundances rather than to solar photospheric abundances, and 
this is the result of the depletion of high-FIP elements rather 
than the enrichment of low-FIP elements. When comparing their 
coronal abundances to solar photospheric abundances, the FIP 
effect persists, but now would be due to the enrichment of low- 
FIP elements, while high-FIP elements remain essentially un¬ 
changed. 
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Parameter: Fe 


Fig. 5: Example contours of overall normalization representing 
EM and Fe abundance (metallicity, as in Fig. [3]> for i Hor (top) 
and HR 7291 (bottom). Abundances were constrained much 
more by averaging across multiple observations. 


5.2. FIP bias and dependence on spectral type 

Continuing our discussion in a more quantitative manner, one 
can define a FIP bias, for instance, inlNordon & Behar|(|2008)l, 


FB = log 


< ^low > 

< ^high > 


(1) 


where < A i ow /hi g h > are the average abundances of low- (< 
10 eV) / high- (> 10 eV) FIP elements. Hence, positive/negative 
FB values reflect a FIP/IFIP bias. FIP biases are summarized 
in Table [6] Errors are calculated assuming Gaussian distribution 
and symmetrizing the errors (average of errors if not originally 
symmetric). We provide three measures here - relative to solar 
values (using all coronal measurements), relative to existing pho- 
tospheric values (using only elements measured in both corona 
and photosphere) and a measure used by |Wood & Linsky ( |2010] > 
(with a negative sign) restricting the FB to the ratio of C, N, O, 
and Ne (or the available subset) over Fe (solar reference). 11 LMi 
and r Boo are consistent with FB - 0 in all measures. On the 
other hand, i Hor , HR7291, and a Cen A and B exhibit signifi¬ 
cant FIP bias with FB > 0 consistent with Fig. [7] The FIP bias 



Spectral Type 


Fig. 6: FIP bias values for the current sample plotted against 
spectral type, using only C, N, O, Ne, Fe as defined in |Wood| 
|& Linsky| ( |2010| ), only negative. Objects marked with an X were 
already included in Wood & Laming (2013). Solid lines bound 
all of the stars in Wood et al. ( 201 2[i except for one. 


with respect to photospheric abundances introduces more uncer¬ 
tainties coming from the photospheric uncertainties, while the 
other two measures generally have smaller uncertainties because 
no error on the solar abundances is assumed. Photospheric mea¬ 
surements with only one reference suffer from the same short¬ 
coming. 

The present FIP biases from Col. 3 in Table |6] are compared 
with the trend reported by [Wood & Linskyj ( ]2010j ) in Fig. [6] Note 
the definition here is the negative of the one defined there. In that 
paper a striking correlation of FIP bias with spectral type (Fig. 9 
therein) is observed for main-sequence G, K, M dwarfs, which 
was later extended in [Wood et~aLl ( |2012| l and [Wood & Laming 
( |2013 i. In the present paper spectral types fall within the range 
studied in those papers, which have already included r Boo , and 
a Cen A and B. Though no obvious trend is apparent for the 
present sample, five of six (at least marginally) fall within the 
FIP bias-spectral type correlation band (see Fig. [6]). The outlier 
is r Boo (spectral type F6), as already noted in Wood & Laming 
P0I3[ along with Procyon (F5), which lies even farther away 
(FB = -0.1). On the other hand, HR 7291, which is of simi¬ 
lar spectral type as Procyon and r Boo (F8), falls nicely on the 
correlation (Fig. [6]). 


6. Conclusions 

We compared the coronal abundances of six stars with their su¬ 
persolar photospheric abundances. For three of them, 11 LMi , 
i Hor, and HR7291, we measured absolute coronal abundances 
from their X-ray spectra. 

- Coronal abundances of all stars are depleted compared to 
their respective photospheres. Since photospheric abun¬ 
dances are significantly higher than solar, depletion is much 
greater than when comparing to solar abundance values. 

- Four of six stars show a positive FIP bias, while the other 
two have no FIP bias. When a FIP effect is present, it is 
different from the solar FIP effect. In the present sample it 
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appears that the high-FIP elements are depleted, while in 
the solar corona it is generally accepted that the low-FIP 
elements are enriched (but this is still disputed, e.g., |Schmelz 
let al.|2012) . 


- Five of the stars (all but r Boo) are consistent with a 
correlation of FIP bias with spectral type (jWood & Linsky| 
20101 [Wood & Laming|20f3) . 


- The importance of knowing the actual photospheric abun¬ 
dances was exemplified by a Cen A and B. In these two, the 
FIP bias is the same regardless of photospheric reference. So¬ 
lar assumptions misleadingly imply enrichment of low-FIP 
elements, while in fact the high-FIP elements are depleted. 
This would lead to different fractionation trends despite the 
similarity in FIP bias. 


In conclusion, high-fidelity measurements of photospheric 
abundances for additional stars with well-determined coronal 
abundances would be highly desirable. Coronal abundance mod¬ 
els should try to explain the phenomena of FlP-independent 
overall coronal metal depletion in stars with supersolar abun¬ 
dances, as well as a more significant high FIP element depletion. 
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FIP[eV] 

Fig. 7: Elemental abundances in solar units of the six coronae (black symbols) as a function of the first ionization potential (FIP) compared with the photospheric values (red). 
Line at 1 represents solar (Asplund et a l. 12009) . Note photospheric abundance ranges represent literature ranges and not actual uncertainties. A FIP effect is observed for i Hor, 
HR 7291, a Cen A and B, which is not true for 11 LMi and r Boo. Clearly, abundance trends relative to solar and relative to actual photospheric values are not always consistent. 
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Table 6: FIP bias measures (uncertainties) 


Object 

Photospheric 

Reference 


Solar 

Wood & Linsky 

O' 

O 

<N 

11 LMi 

0.08(0.38) 

0.00(0.43) 

0.08(0.44) 

l Hor 

0.38(0.20) 

0.31(0.18) 

0.54(0.18) 

HR 7291 

0.50(0.41) 

0.32(0.17) 

0.92(0.17) 

r Boo 

0.42(0.75) 

0.23(0.54) 

0.28(0.38) 

a Cen A 

0.45(0.22) 

0.37(0.17) 

0.65(0.16) 

a Cen B 

0.34(0.20) 

0.28(0.16) 

0.42(0.18) 


1 A-Fel < Ac,Ao,An,Ani; > 
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